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Abstract 

We present a systematic determination of the order h fermionic energy shift 
when an electroweak string is perturbed. We show that the combined effect of zero 
modes, bound states and continuum states is to lower the total fermionic ground 
state energy of the string when the Higgs instability of the string is excited. The 
effect of the Dirac sea is thus to destabilise the string. However, this effect can 
be offset by populating positive energy states. Fermions enhance the stability of 
an electroweak string with sufficiently populated fermionic bound states. 
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1. Introduction 

The resurrection of interest in electroweak strings, first discovered by Nambu 
[0, induced by the discovery of the W and Z-string solutions lead to intense 
research into the properties and potential impact of these objects. Both the W 
and Z-strings solutions are U(l) Nielsen-Olesen vortices embedded in some U(l) 
subgroup of the Weinberg-Salam model [Q]. Unlike the archetypal Abelian Higgs 
Nielsen-Olesen vortex and standard GUT vortices, the stability of these embed- 
ded strings is not guaranteed by any topological argument. The enlarged field 
space of the SU(2)xU(l) Salam Weinberg model provides extra directions in field 
space which allow the construction of continuous, finite energy paths that unwind 
the string solution. The criterion for the existence of topological strings is that 
the first homotopy group of the vacuum manifold is non-trivial [^. The first 
homotopy group of the three sphere, which is the vacuum manifold of the elec- 
troweak theory, is however trivial and therefore the electroweak string solutions 
are non-topological in nature. This implies that the W and Z-strings are just 
finite energy solutions to the classical field equations and, at best, electroweak 
stings could sit at a local minimum of the electroweak energy functional. They 
would then be separated from, and prevented from decaying into, the trivial vac- 
uum by a finite energy barrier. This barrier would either have to be surmounted 
or quantum mechanically tunnelled through to destroy the string solution. Alter- 
natively electroweak strings could be local maxima or saddle points of the energy 
functional, in which case the pure string solutions would be intrinsically unstable. 
The issue of stabihty is a question of energetics and depends on the value of the 
parameters used. Detailed stability analyses of the Z-string have been carried 
out by two independent sets of authors[^, |^. Both investigations concluded that 
the Z-string solution possess modes of instability except in the unphysical region 
where siri^i^Ow) ~ 1- Other investigations imply that the W-string solutions 
are unstable for all values of the electroweak parameters. W-strings will not be 
considered any further during the course of this paper. 

The instability of electroweak strings for physical values of the Standard Model 
parameters would tend to imply that electroweak strings are physically insignif- 
icant. Had electroweak strings been stable then we would have expected a net- 
work of strings to be formed during the electroweak phase transition through the 
Kibble mechanism. The density of electroweak strings would be lower than their 
topological counterparts as the vanishing of the electroweak Higgs field at a point 
does not necessarily imply the existence of a string. Another important differ- 
ence arises because electroweak strings would tend to form fragments of finite 
length that are bounded by a monopole anti-monopole pair. As discussed in [|| 
the inherent contraction of these string fragments within the electroweak theory 
provide all of the conditions required for electroweak baryogensis. Although the 
CP violation in the minimal electroweak theory is too small to account for the 
observed baryon anti-baryon asymmetry, electroweak string are also expected to 
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arise in extensions of the minimal electroweak which admit greater CP violation 

i- 

The stability analysis discussed thus far has focussed on the pure Z-string solu- 
tions. It is however possible that the formation of bound states of other particles 
on the string may have a non-trivial effect on the stability of electroweak strings. 
Of particular importance for a string like solution is the effect of a finite density 
of bound states moving along the string. In |jlO[| the authors considered a toy 
model in which they extended the electroweak theory by adding an extra global 
U(l) scalar to the bosonic sector of the electroweak theory. They constructed 
solutions with a finite U(l) charge per unit length of the string and then investi- 
gated the stability of the resulting solutions. They found that it was possible to 
extend the region in which the electroweak string is stable to si'n?{9w) ~ 0.47. 
Thus the stability of an electroweak string may be enhanced by its interaction 
with other fields. The stabilising effects of a plasma of scalar particles have also 



been considered \\\ 



The authors of |[TD| rightly point out that, in the context of the electroweak 
theory, the pertinent particles to consider are the fermions which gain their mass 
via a Yukawa coupling to the Higgs field. As will be discussed in this paper, 
the systematic analysis of fermionic states in the electroweak string background 
is a far more involved task than the analysis of bosonic condensates. The crit- 
ical difference arises because of the different statistics satisfied by fermions and 
bosons. Bose-Eienstein statistics allow a significant population of the bosonic 
ground state and hence a classical treatment of this single mode is appropriate. 
In sharp contrast, the Pauli exclusion principle does not allow a similiar clas- 
sical treatment in the fermionic case and a complete field theoretic analysis is 
required. Naively, as the fermions gain their mass through a Yukawa coupling to 
the Higgs field it would be expected that the massive fermions would rather be 
localised about the centre of the string and prevent the string from decaying in 
an attempt to maintain their low energy state. The purpose of this paper is to 
carryout a systematic analysis of fermions in the background of the Z-string. The 
effects of the fermionic zero mode on the stability of the string was considered 
in |T2| , pTSf where it is claimed that the fermionic ground state energy makes a 
negative contribution to the energy of the string and increases the instability. 
However, we must consider the full Dirac spectrum. It is not clear what effect 
the infinite number of continuum states have upon string stability and positive 
energy bound states can only be consistently included if the actual numerical 
value of the fermionic ground state energy is calculated. Here we present a sys- 
tematic, quantitaive treatment of the fermionic contribution to the energy of the 
Z-string. 

It will be shown that a systematic treatment of fermions in a string back- 
ground requires a complete field theoretic calculation. Techniques are developed 
for the computation and renormalisation of the Dirac sea contribution to the 
energy of the Z-string. Many of the techniques developed have a wider range 
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of applicability, particularly for determining Dirac sea energies in other 3+ ID 
backgrounds such as the sphalcron. Having shown how the Dirac sea energy can 
be computed in the Z-string background, the change in this energy is computed 
when a mode of instability of the Z-string is excited. This determines the con- 
tribution of the fermionic ground state to the stability of the Z-string. Having 
understood the ground state energy it is possible to consistently consider the 
effects of filled, positive energy bound states on the stability of the string. 

In section 2 we review the electroweak string and its stability. A systematic 
discussion of the semi-classical treatment of fermion-soliton interactions is given 
in section 3. In this section the origins of the Dirac sea contribution to the 
fermionic energy of a soliton are demonstrated from a field theoretic point of 
view and a formal statement of the required calculation is made. In section 4 we 
review the Dirac equation in the background of the Z-string and derive the whole 
fermionic spectrum. 

The fermionic ground state energy is divergent and its computation requires 
the use of a regulator and renormalisation scheme. In section 5 the proper time 
regularisation procedure is applied and an expression for the regulated fermionic 
energy of the pure Z-string is derived. We present a consistent renormalisation 
of the regularised fermionic energy in section 7. The appropriate counter terms 
are derived using the heat kernel expansion. 

In order to calculate the fermionic energy, we must sum over all continuum 
states. This sum is implemented by discretising the continuum, determining mo- 
mentum shift functions and analytically recovering the contimmm. The formal- 
ism for this procedure is established in section 6. To demonstrate the consistency 
of this procedure and highlight its key physical properties, we analyse a core 
model string in section 8. The core model introduces simplified string profile 
functions which allow the Dirac equation to be solved analytically and allows sig- 
nificant insight into numerous subtle aspects of the calculation of the fermionic 
energy to be gained. 

The calculation of real physical interest is presented in section 9. We in- 
troduce the perturbed string, determine the Dirac spectrum in this background 
and evahiate the counter terms. We then apply the formalism of sections 5 to 
7 to determine the difference between the ground state fermionic energy of the 
perturbed string and that of the pure Z-string. 

Once the change in the fermionic ground state energy has been calculated it 
is only a small step to add a finite density of positive energy bound states moving 
along the string. In section 9 we investigate the effect of changing the density of 
positive energy states upon the stability of the string and discuss the possibility 
of stabilising the electroweak string. 

We conclude with a summary of what has been learned during this investiga- 
tion and point towards future avenues of research. 
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2. Strings in the Electroweak Theory 

In this section we review the Z-string solution and its stabihty. The existence 
of string in the electroweak theory was first proposed by Nambu 0] in 1977. 
Interest in electroweak string was however resurrected by Vachaspati with the 
discovery of the Z-string and W-string solutions. 

To set notation, we take the purely bosonic sector of the electroweak theory 
to be defined by the Lagrangian density, 

LBosonic = --^W^^W^^" - \f^.F^^''+ I |2 -A(| $ p -ri^f (1) 

where, 

= d,W: - d.W; + ge^^'Wlw:, (2) 

d^B, - d,B^, (3) 



D,^ = {d, - - 1-5,)$, (4) 

and, 

is the standard Higgs doublet. The gauge boson mass eigenstates are defined in 
the usual way, 

(l)-{:zl: -it)!":)' -^->.'±-.^). 

where Oyj is the Weinberg angle, g = 2gcos^^ and g' = 2qsin6.^. 
The Z-string ansatz is defined by 0, 

^^ = _^sin^, Z2 = ^cos^, <p^ = r]f{r)e'\ (7) 
qr qr 

with all other fields in the theory set to zero. By keeping only the fields which 
are non-zero in the Z-string background, the energy of the field configuration 
becomes, 

E = j (fx {^Zl+ I d,(t>o - ^qZ4o r +A(| 0„ I' , (8) 

where the indices i, j = 1, 2, 3. Rescaling the radial variable by r — > M^r, allows 
the energy per unit length of the string (in units of M|) to be written as, 

E.tr^n, = 2n f dw ( (^) ^ + {ff + [^) ^ " 1)' + - , (9) 



where, (3 = Variation of (|^) with respect to the profile functions /(r) and 
i/(r), leads to the familiar Nielsen-Olesen equations for a unit winding number 
string, 

r (r) + ^(Kr) - 1)' - Pifirf - I) fir) = 0, (10) 



r 



^"^r)-'^-{u{r)-l)f{r)=0. (11) 
r 

Thus the Z-string is a U{1) Nielsen-Olesen string embedded in the electroweak 
theory in the direction of the U{l)z subgroup generated by the Z gauge particle. 
The Z-string is not topologically stable, such a field configuration can decay 
to the trivial vacuum, and stability is an issue of energetics. By searching the 
(/?, sin^^^) plane it is found that the electroweak string in unstable except for 
sin^ 6'^ ^1 ^. Of particular importance, electroweak strings are unstable for 
the physically observed values of the electroweak parameters ( sin^^^ ~ 0.23). 

For a unit winding number string two modes of instability were found. The 
first corresponding to the upper component of the Higgs doublet, 0+, acquiring 
a non-zero value in the string core and the second corresponding to the W field, 

pie j 

= -^{W- + -Wi), (12) 

forming a non-zero value inside the string core. The latter instability corresponds 
to the formation of a W-condensate inside the string, this mode of instability has 
been investigated in [TB|. The effect of these modes is to cause the string to 



unwind and decay to the trivial vacuum. 

The stability of electroweak string may however be changed by the presence of 



bound states. Vachaspati and Watkins found that, by coupling the bosonic 



sector of the electroweak theory to a global complex scalar field, it is possible 
to significantly enhance the stability of the Z-string solution by the formation 
of scalar bound states. However, as they point out, the pertinent particles to 
consider are the fermions of the standard model. As the Higgs field vanishes in 
the core of the string, fermions which gain their mass from a Yukawa coupling to 
the Higgs field will be massless inside the string core. Thus naively they would 
rather reside within the core of the string and resist its dissolution. 

While the scalar case can be considered classically, a purely classical analysis 
is not possible for fermions as their statistics prevent a large occupation number 
within a given energy state. Consequently a full quantum mechanical treatment 
is required. 

Fermions in the background of the electroweak string have been considered 



by several authors |T^, |T2|, |T^]. However, only the effects of the zero mode states 
were considered. We present a systematic calculation of the change in the entire 
fermionic energy brought about by perturbing the Z-string. 
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3. Fermion-Soliton Interactions and the Dirac Sea 



The fermionic part of the electroweak action can be written in the form, 

i J ci^x^U[$,Z^]^ (13) 

where Z^] is some operator which depends upon the background scalar and 
gauge fields and the fermionic fields are grouped together in the spinor \E'. The 
fermionic contribution to the path integral is a Gaussian functional integral over 
the fermionic fields. For anti-commuting fields this Gaussian functional integral 



has the standard result ITB 



J D[^^{x, t)]D[-^{x, t)] exp (^i J Z^]^^ = Det i[$, Z^]. (14) 



Using (|Tj) it is possible to formally integrate out the fermionic contribution to 
the action and obtain an effective action which encapsulates the effects of the 
fermion fields yet only contains explicit dependence on the bosonic fields. This 
can rarely be done in practice because the calculation of the determinant in (0) 
requires a knowledge of all the eigenvalues of the operator Z^] for a general 
background. Instead we consider specific backgrounds. 

The operator Z^] is related to the Dirac Hamiltonian Z^] in the 
bosonic background of the fields $ and Z^ by, 



A[^,Z,]=z^^-H[^,Z,]. (15) 



For a static field configuration the Dirac Hamiltonian has the set of energy eigen- 
values {Er} which satisfy, 

H[<!>,Z,]^r = Er^r. (16) 

This is just the Dirac equation describing the fermions in the background of the 
string. 

The determinant in ([1^) can be expressed in terms of the energy eigenvalues 



IS 



iT (j2i-Er + nrEr 



(17) 



Det A[$,Z^]= ^C(K})exp 

{rir} 

where the integers are the occupation numbers of the excited positive energy 
states and C{{nr}) are the associated combinatoric degeneracy factors. 
The fermionic ground state energy is then given by, 

E,s = Ei-Er). (18) 
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This is exactly the energy which would arise from filling all the negative energy 
states as predicted by Dirac. 

The ground state energy of the field configuration is then given by, 



E ground EQinggi^iii -\- ^^Eiyos ~\~ ^ ^ ( ^r)^ 



(19) 



where E^os is the effect of the bosonic fiuctuations (these will not be considered 
here) . 

We will be interested in the difference between the Dirac sea energies of two 
systems, this will generally be referred to as the fermionic energy difference, which 
can be written as, 

SEfermion = ("^'■.l ~ -^^-.a) ) (20) 
r<0 

where and Er^2 are the energy eigenvalues of the Dirac operator in the two 
backgrounds and the sum is over all the negative energy states. It is useful to 
express the fermionic energy difference ( pO|) as a functional trace over a complete 
set of eigenf unctions. Specifically, 

= -^Tr (I i^i I - I ^2 I) = -^Tr (^/^- , (21) 

where Hi^2 are the Dirac Hamiltonians for the two systems and the factor of | 
arises because the trace in (pTf) is over both the positive and negative energy 



eigenstates. These are known to be in a 1 : 1 correspondence through the action 
of the charge conjugation operator. 

The computation of ( pTD is the aim of this paper. All energies per unit length 
of string will be expressed in terms of the square of the Higgs VEV. 

4. The Dirac Equation in an Electroweak String Background 

Fermions enter the electroweak theory through the Lagrangian density, 

L fermion = I^lI^D^^Il + IUr'^^D^Ur + idRj^D^dn + Lyukawa, (22) 

where, 

Lyukawa = -9u(^L^Ur + M/J^^^l) " 9d(^L^dR + rf/j^^^i), (23) 

is the Yukawa term which is responsible for generating the masses of the fermions 
after symmetry breaking. 

$ is the standard Higgs doublet defined by (^ and $ is defined to be. 



The left handed fermions are combined in the usual SU{2) doublet defined 

by, 



Ul 
dL 



The covariant derivatives, which couple the fermions to the gauge fields, are 
defined by, 

D^^L = {d, - i^-a'^W; - i^YLB^)^L, D^^r = {d, - 

where Yl and Yr are the f/(l)y hypercharges and g, g are the SU{2) and U{1)y 
coupling constants respectively. Throughout explicit reference is made to the up 
and down quarks through the notation ul,r and dL,R- Other fermions can be 
included by simply introducing the appropriate hypercharges. 

Varying the Lagrangian with respect to the fermion fields produces four two 
component Dirac equations, 

h'^d^UR - qaR-i^Z^UR + ee^'^^A^^UR - gu(f>oUL + gu(t>+dL = 0, 
il^d^^UL + qal-f^Z^UL + ee''-f''Af,UL + e^YW~dL - gu(f>*oUR - gd(t>+dR = 0, 
il^d^^dR - qa%rt^Z^dR + ee'^^^^^dR - gd(f)*+UL - gd(t>*odL = 0, 

il^d^UL + qai-ff^Z^dL + ee^^YA^dL + e^jf^W+UL + gu(t>*+UR - gd(t>odR = 0, 

where, e = g sin 26*^4, and e^^ = \/2qcos0y,. 

Throughout we use the Wcyl representation of the Dirac matrices. 
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-1 
-1 



;7 



a' 
-a' 



i7 =75 



1 
1 



with the chiral projection operators defined in the standard way, 
Pr 



^(1 + 7^), ^L = ^(l-7'). 



Working in temporal gauge (W^ — Ag — Zg — ) and isolating the time 
dependence implies that the Dirac equations read. 



idf 



/ Mr \ 


f 


Ul 




dR 




\ dL ) 


\ 







gu(p+ 



\ 


( ur\ 




Ul 




dR 


I 


\ dL J 



(24) 
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where i = 1,2,3 and the covariant derivatives are defined by. 



Dfi? = di + iqa^^Z, - tee^Ai, Df^ = di - tqa^Zi - iee^Ai. 

In the background of the pure Z-string the only non-zero bosonic fields are Zi 
(where i = 1,2) and 0o, the specific forms of which are given in (^. The Dirac 
equation (^) can be seen to decouple into two independent four component 
equations for the up and down quarks. 

As the Z-string solution is independent of z, the Dirac equation in the back- 
ground of the Z-string can be written as, 

tdtP - iCd^P = H2+i,xP, 
where P denotes a general 4 component spinor and. 



H2+i,\ is the Dirac Hamiltonian defined in the plane perpendicular to the direc- 
tion of the string, i takes values 1, 2, and the parameter A allows much of the 
analysis for the up and down quarks to be carried out using the general Hamil- 
tonian (|25D , with A = 1 for the up quark and —1 for the down quark. The U{l)z 
gauge charges then satisfy, 

+ ftR = A. (26) 

The parameter A arises as the left handed up quark couples to (po and the 
right handed up quark couples to the complex conjugate field 0*, whereas the 
roles are reversed for the down quark with the left handed component coupling 
to 0* and the right handed component coupling to 0o- 

The full 3 + ID Dirac Hamiltonian can be written as, -ff3+i,A = H2+i,\ + 
iCS^. Setting, 

P(x,t) = P{xi,X2)e'^'+''e-'^'\ 

where -Es+i is the total fermionic energy and kz is the momentum along the 
direction of the string, the square of the Dirac equation becomes, 

^3+l,A^(^l'^2) = + klP{xi,X2), 

where we have used {C, H2+i,\} = 0. The problem thus reduces to solving the 
2 + ID Dirac equation, 

E2+i,xP{xi,X2) = H2+l,xPiXi,X2). (27) 

Using {C, H2+i,x} = 0, if there is a solution of the 2+1 dimensional Dirac 
equation, P{xi,X2), with energy -E2+i,a > then CP{xi,X2) is easily seen to 



10 



be a solution with energy — £^2+i,a- The particle conjugation operator therefore 
provides a 1 : 1 mapping between the positive and negative energy solutions of 
the 2 + ID Dirac equation. 

Introducing polar coordinates in the plane perpendicular to the direction of 
the string, the 2 + ID Hamiltonian becomes, 



2+1, A 



—ze 



ie 







V 

















le 



-ixe 



-le 





dr-idg 





which acts on a spinor of the form. 



V J 



E 



The radial profiles of the spinor components are then determined by. 



{E3+i,\ + K)^)^ 



dr - - + aR—^ j - 9f{r)ipA, 



dr 



h Ql ^4 - gf{r)ti, 



(28) 



(29) 



(30) 
(31) 
(32) 
(33) 



At small r, the Nielsen Olesen profiles give /(r)~r and i/(r)~r^, which give 
the following small r forms for the spinors. 



/ 1 \ 




V / 



1 


V / 



,(n-A) 



/ \ 



1 

V J 



-(n-A+l) 







V 1 / 

(34) 

At large r, the Nielsen Olesen profiles give /(r)~l and i/(r)~l, which give 
the following large r forms for the — equations of motion. 



Eipi = -(dr + ^ ) i)2 - gips, 



(35) 



11 



-#4, 

V'4 - #1, 



where, 



^ n - X + aL- 



(36) 
(37) 
(38) 

(39) 



The spectrum of states consists of a continuum, bound states and zero modes. 

If E'^ > g"^, the four hncarly independent regular solutions at large r can be 
expressed in terms of the Bessel functions Ju{z) and N^{z), 



( '^1 \ 

^3 
V ^4 / 



+03 





V ±fJ|.,|±i(A;r) J 

( (kr) \ 

+1 A^|..|±i (kr) 




I 



+ 02 







fe"^I^Al±l 

{kr) 



[kr) 



\ 



a4 



V ±f ^|..|±i {kr) j 

i ^ 

{kr) 

{kr) 

V ±f {kr) J 



(40) 



where, E'^ — k"^ -\- g"^ , is the usual mass-shell condition. The upper signs corre- 
spond to the case when \ i'\\= fx and the lower signs when | z^a |= —I'x- 

If E"^ < g"^, the asymptotic radial solutions are modified Bessel functions and 
the general solution can be written as. 



V ^4 y 



/ {kr) 



ai 



■f ^-A+i (kr) 




+ 02 



+03 



V ilu,+i{kr) J 

( K,. {kr) \ 
{kr) 



V -f {kr) J 





-f {kr) 
{kr) 



\ 



V V^,+i{kr) J 
( 



+ 04 



(fcr) 



(41) 



where, E"^ ^ g"^ - /c^, with /c^ e (0,^^). 

Only the last two independent solutions are square integrable at large r, thus 
a physical bound state solution must match onto these solutions only. 

For each value of n there arc two linearly independent, regular short distance 
spinors, each one matches onto a unique linear combination of the four asymptotic 
large distance spinors. We denote the matching coefficients by as,u where s = 1, 2 
labels the short distance spinor and / = 1 — 4 labels the asymptotic spinor. To 
produce a normalisable state a suitable linear combination of the two regular 
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0. (43) 



short distance spinors must be taken such that the coefficients of the divergent 
asymptotic solutions are zero. Taking the hnear combination, aSi + (3S2, of the 
two regular short distance spinors, square integrabiltiy implies, 

ai^ia + a2,i/5 = 

01,20 + 02,2/5 = 0. (42) 

A non-trivial solution to (|42|) exists if and only if the associated determinant is 
zero, that is, 

' ai,i (32,1 

(^1,2 0,2,2 

In general, finding bound state solutions is reduced to the problem of finding the 
values of the energy which are roots of (^). 

A special case are zero modes |T^. Setting E = in the equations of motion, 
we see that the system decouples into two second order systems for all r. 

At small r the leading order behaviour of the two solutions to the first set is 
given by, 

^i^r", ^4~r("+i) and ~ r'^""^), ^4 ~ 
While the second system gives, 

^2 ~ r~("+^), V's ~ r"", and 7/^2 ~ r("-^+^), 7/^3 ~ r^""^). 

Both systems possess one regular and one irregular asymptotic solution. Thus 
if both short distance solutions in one system are regular, they will match onto 
one regular and one irregular asymptotic solution. In the language used above, 
cti,2 = CL2,2 = (say), one constraint is automatically satisfied and there is always 
a linear combination of short distance spinors that yields a square integrable 
solution. Thus we have zero mode solutions. 

By setting A = 1 it is easily seen that the only case in which there are two 
regular short distance solutions in one system is for n = 0. The Up quark zero 
mode therefore has n = and V'2 = V^s = 0. Further we see we still have a 
solution to the equations of motion with '?/'2 = V's = if we set -Es+i.a = ^z- 
Physically this means that the up quark zero mode can move along the positive 
direction of the string at the speed of light. This should be compared with the 
massive bound states which are free to move in either direction along the string. 
This restriction on the motion of the zero mode has important implications for 
the formulation of the Dirac sea energy. 

Similarly, by setting A = — 1 we find the only down quark zero mode if n = — 1. 
The solution has ipi = ip^ = and -E3+i,a = —kz, thus this mode is restricted to 
move in the negative z— direction at the speed of light. 

The existence of these zero mode solutions in the pure Z-string background is 
guaranteed by an index theorem PD[. If the string is perturbed, the zero modes 
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are lifted and become (for small perturbations) low energy massive bound states 



1^, |T3|- The existence of massive bound state solutions is not guaranteed by any 



form of index theorem but depends upon the precise values of the parameters in 
the theory. For each set of parameters, solutions to the bound state condition 
l3f) must be searched for. 



5. Regulating The Fermionic Energy 

We have seen that the fermionic vacuum energy can be calculated by summing 
over all the negative energy eigenvalues of the Dirac equation in the background 
of interest. Clearly this quantity is divergent: there are an infinite number of 
continuum states which make progressively larger contributions to the sum. This 
problem is closely linked with the usual ultraviolet divergences which arise in 
quantum field theory when large values of loop momenta are integrated over. 
Effectively the Dirac sea energy corresponds to a 1-loop calculation in the back- 
ground of some field configuration. As with normal loop calculations these di- 
vergences can be consistently handled using a renormalisation procedure. The 
issue of regulation is discussed in this section and renormalisation is discussed in 
section 7. 

It is useful to look in a semi- quantitative manner at the nature of the ultra- 
violet divergences expected. The continuum energy contribution takes the form 
of an intergal over all the momentum components, the nature of the divergences 
can be isolated by introducing a naive ultraviolet cut-off. A, 

E,ont ~ j'^d^kE ~ 0{A^) + 0(A2) + O(logA) + ■ ■ ■ 

The first term above will cancel when two strings are compared, leaving quadratic 
(A^) and logarithmic (log A) divergences which must be dealt with through the 
renormalisation procedure. 

The freedom to move along the length of the string also means that the bound 
states contribute to the divergences in the fermionic energy, 

Etound ~ J^'dk^Mi + k^ ~ 0(A2) + 0(log A) + • • • (44) 

When two string are compared it is found that there is a natural way of pair- 
ing the states of the two systems. The result of pairing the states is that the 
quadratic divergence in (|^) cancels when two system are compared and the re- 
sulting divergence is logarithmic. It is this logarithmic divergences which must 
be removed by the renormalisation procedure. 

These divergences are of exactly the types expected, as standard Feynman 
diagram calculations in the electroweak theory produce both quadratic and log- 
arithmic divergences. The ultraviolet behaviour of a theory is associated with 
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very large momentum and thus probes very small distances hence the types of 
the divergences expected should be independent of any external field configura- 
tion. 

As discussed in section 3, we are interested in the functional trace, 

6Ef,rmion = "^Tr (^/^ - , (45) 

where the two backgrounds being compared have Dirac Hamiltonian operators 
Hi and H2- To regulate the divergences we introduce a proper time regulator, 

5^/e™(r) = r ^ Tr (e-*^? - e"*^^'). (46) 

4\/7r Jt t2 

T acts as an ultraviolet cut off through the exponential suppression of high mo- 
memtum contributions and we are interested in the t —>■ limit. It can be seen 
on dimensional grounds that the proper time regulator, r, is related to the naive 
ultraviolet cut off. A, by, r ~ 

We can see that this regulated form reproduces the original trace in the limit 
r ^ by using the standard integral pT| , 

r =1 ^ I E'r) -^-2V^\E I +0(r^). (47) 

It is important to note that two systems must be compared so that the ^ term 
in (^) cancels. 

The spectrum of the Dirac Hamiltonian in the background of the Z-string 
consists of three distinct parts: zero modes, bound states and continuum states, 
thus the regularised fermionic energy can be written as, 

and each of these three parts can be handled separately. 

As discussed in section 4, we split the 3 + ID Dirac Hamiltonian into a 2 + ID 
Hamiltonian together with a contribution coming from the momentum of the 
fermion along length of the string, -fff+i = -ffl+i + ^z- The trace must sum over 
all the allowed values of kz, remembering that zero modes only move in one 
direction along the string. We restrict the spinors to a cylindrical box of length L 
and impose periodic boundary conditions in the direction, thus quantising the 
allowed ^-momenta. Taking the limit L oo, the sum over allowed momenta 
becomes an integral over the momentum parallel to the string in the usual way. 
As a zero mode only moves in one direction along the string, its contribution is 
only half that of a massive bound state or a continuum state. 

We can now determine the zero mode contribution to the energy per unit 
length of the string; 

^ ^_ /-^ rft p rffc _4fc2 _ OzM /-^ ^ _ OZM 1 
^zeromode\T) /-^ n /— r, ^ o /j-2 o ' \ °) 
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where azM denotes the number of zero mode solutions to the 2 + ID Dirac 
equation. From (|48|) it can be seen that a zero mode state contributes a purely 
manifest quadratic divergence to the fermionic energy. 

Similarly for a massive bound state solution with effective mass E, the energy 
contribution per unit length of string is given by, 



where C is Eulers constant. As with the zero modes it can be seen from (^) 
that each massive bound state contributes a manifest quadratic divergence to 
the fermionic energy. Moreover in both cases the manifest quadratic divergence 
is independent of the energy of the state. We shall see in section 8 that these 
manifest quadratic divergences cancel when two systems are compared as there 
is a natural way in which to pair states. 

Finally the continuum contribution to the fermionic energy per unit length of 
string reads, 



where Tr^„j denotes a trace over the positive energy part of the continuum spec- 
trum of H2+1 and we have integrated out the degree of freedom along the direction 
of the string. 

The contributions to Efermionij) from the massive bound states and zero 
modes are easy to compute once the fermion spectrum is known. The continuum 
contribution in contrast requires a significant amount of work before it is of a 
practical form for calculations. 

6. Discretising and Recovering the Continuum 

In the previous section we determined an expression for the regulated contri- 
bution of the continuum to the total fermionic energy in the string background. 
This involves a trace over the positive energy continuum and consequently in- 
troduces a sum over an infinite number of states. To carry out this sum we 
impose a boundary condition on the spinor profile functions on the surface of a 
cylinder of radius Ro centred on the core of the string. This condition selects a 
discrete subset of the continuum states over which we trace. The full continuum 
is recovered as i?o 00, so the actual trace over the continuum can then be 
generated by taking this limit. This will be done explicitly so that the discrete 
sum over the discretised continuum becomes an integral. The use of the dis- 
cretisation condition is merely an intermediate step in setting up the sum over 




(49) 




(50) 
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the continuum. This step is however critical, because it ensures that the correct 
density of continuum states is used. 

The details of the discretisation procedure are given in appendix A. The 
boundary condition we impose is that the 1st and 3rd components of the spinors 
vanish at r = Rq. 

The Dirac equation in the Z-string background has two regular short distant 
solutions which match asymptotically onto two linear combinations of the asymp- 
totic solutions given in (0). The boundary conditions lead to a constraint of the 
form, 



J\u\{kRo) 



-X{k), (51) 



where X{k) is given in terms of the matching coefficients. 

The left hand side of this constraint is approximately cot(fc-Ro) and consists 
of an infinite number of branches of width proportional to The density of the 
solutions to the discretisation condition thus becomes infinite as -Ro ^ oo and so 
the full continuum is recovered. 

Obtaining a unique solution to the discretisation condition requires careful 
analysis, principally due to the occurence of singularities in X{k). This issue is 
discussed in detail in appendix A. 

The discretisation condition can be written uniquely in the form, 

hRo = k^Ro + ^iki), (52) 

where the free momentum is defined by, k°Ro = iir, and A{k) is determined by 
matching short distance spinors to asymptotic spinors. The details are given in 
appendix A. 

As, ki = k° + 0(-^), we have, 

kiRo = KRo + A±(fc°) + O(^). (53) 

All of the terms on the right hand side depend only on the known free momentum 
k°. Expanding the momentum in powers of we have. 



-tk^ -* 
e ' = e 



2k° 1 
Ko it 



(54) 



The free momentum k° provides a set of evenly spaced nodes with spacing 
6k = In the large Rq limit we have, 

V e-*^' = ^ / dke-"' - ^-fl dkkA^{k)e-"'" + 0(-^). (55) 

~^ IT Jo 2 TC Jo Ro 

The first term on the right hand side of (^) will be seen to cancel when two 
systems are compared. 
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Equation |55| is exactly the integral representation of the continuum contribu- 
tion to the fermionic energy we require. We have expessed the regularised trace 
over the continuum states as an integral over the transverse momentum with the 
integrand defined in terms of matching coefficients found by solving the Dirac 
equation. 

Finally, the total fermionic contribution to the energy of the string from a 
given angular momentum mode, ra, is given by, 



1 1 dt 



E%rrmon{r) - + j ^ ^2 E ^ 



?2 



27r 





'Ti r 


Jr 


2~ Jo 



(56) 



where cx'^m is the number of zero modes with mode number n, {-En,i} is the set 
of positive energy massive bound state solutions of the 2 + ID Dirac equation 
and finally A^'"'(/c) are the two functions defined by ( |152|) for the pure Z-string 
background. 



7. Renormalisation Issues 

Now we have a regulated expression for the fermionic energy shift we must 
renormalise it in order to calculate the physical value of the fermionic energy. We 
focus on the computation of the counter term in the background of the Z-string 
using the heat kernel expansion. By calculating the counter term required to 
renomalise the contribution made by each angular momentum mode it is possible 
to renomalise each mode individually. 

The highest order divergence which is expected in the electroweak theory 
is a quadratic divergence associated with the Higgs mass renormalisation. The 
remaining divergences are logarithmic in nature, corresponding to the wavefunc- 
tion renormalisation of the Higgs and gauge fields together with the logarithmic 
renormalisation of the Higgs self coupling. 

The physical energy of the string should be expressed in terms of physical 
parameters defined at the electroweak scale, rather than the bare parameters 
that enter the Lagrangian. Consider the sum of the classical string energy and 
the fermionic energy, 

Estring = E^i^gg^^^i + E f (^j-mionij) ) (^7) 

where E^^f^^^^^^i is the classical energy per unit length of the string expressed in 
terms of bare parameters (see (^) and E fermion{T) is the regularised fermionic 
energy. The divergences arising in the fermionic energy should be combined with 
the classical energy defined in terms of the bare parameters to give an expression 
for the classical energy defined in terms of the physical parameters defined at 
some scale /i. 
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In the proper time regularisation scheme we have, 



The integral over the proper time t can be spht into two pieces, t ^ fi and 
yU oo, where /i is some arbitrary scale (eventually /i will be set equal to M^^). 
We use the heat kernel expansion to isolate the divergent contributions arising in 
the small r limit of the first integral above and determine the counter term, 

1 /-/^ dt^ 



EcT{r,fi) = — T^Tre-*^' U, . (58) 



IT Jt t 

Consistency of the renormalisation scheme requires that no physical quantity 
depends on the scale /x. This will be ensured by renormalisation group invariance, 
provided the terms which comprise EcT{T,fJ') take the same form as those which 
arise in -Eclossica/ there exist unique expressions relating the bare parameters 
to the renormalised parameters. These conditions are satisifed in the electroweak 
theory when the proper time regulator is used |Q. Egtring is thus invariant under 



changes in fi and so the value of n will be set to an energy scale characteristic 
of the electroweak theory and the parameters of the theory set to their values at 
this energy scale. Specifically, fi = M^^. 

In the electroweak theory anomolies cancel between the quark and lepton 
sectors. As we consider only quarks, there is the possibility of anomalous diver- 
gences appearing that can only be cancelled by including leptons. However, as 
the two systems being compared are in the same topological sector, no anomalous 
divergences arise. 

To compute the counter term we use the heat kernel method. First we review 
the method and then we apply it in the string background. 
Consider an operator, O, and the associated heat equation, 

^^H{t;x,y) + d,H{t;x,y) = 0. (59) 

The trace of interest is then given by , 

Tr e-*° = Tr y" d'^xH{t; x, y). (60) 

It is important to make the distinction between the two types of traces which 
appear in (pO|), the trace on the left hand side is a functional trace over a complete 
set of states whereas the trace on the right hand side is a standard matrix trace. 

It is well known that the heat equation possess an asymptotic solution in 
terms of the heat kernel expansion defined by p2| , p3| . 



1 I -r _ |2 

H{t- X, y) = exp( l /I ) ^ a^x, y)t\ (61) 



n=0 



19 



( |6TD is a power series expansion for small proper time t. The functions an{x,y) 
are the heat kernel expansion coefficients, their precise form will be derived later. 
Substitution of (|6T|) into ( |60D imphes that, 



Tr e 



-to 



1 



(62) 



(47rt)2 



2 n=0 



where [a„] = a„(x, x). 

By direct substitution of (pl| ) into (|^) and equating powers of t it is straight- 
forward to derive a recurrence relation for the coefficients an- 

The counter term defined at some subtraction point, /i, involves a trace over 
the 2 + ID Dirac Hamiltonian defined in the plane perpendicular to the string. 
Therefore d = 2 in the heat kernel expansion and we have. 



term corresponds to the quadratic counter term and the [02] term correspond to 
the logarithmic counter term. The quartic divergence arising from the [oq] term 
cancels when two systems are compared. 

To derive the counter term only the coefficients [ai] and [02] need to be cal- 
culated. They are found to be. 



where /(r) and z/(r) are the Higgs and gauge field profiles of the string respec- 
tively. 

In this section the forms of both the quadratic and logarithmic counter terms 
have been derived using the heat kernel expansion. The next step is to calculate 
the difference between the fermionic energies of two strings. In the next section 
such a calculation is carried out using a toy model which provides significant 
insight into the calculation of fermionic energies. 




(64) 




(65) 
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8. The Core Model 



Before attempting a full numerical calculation of the fermionic energy in the Z- 
string background, it is very instructive to consider a more analytical treatment 
in a simplified model. To this end we use a core model in which simplified 
profiles replace the actual Nielsen Olesen profiles. Exact expressions for the 
spinor wavefunctions can then be written down, the discretisation condition can 
be constructed analytically and insight into some of its important properties can 
be gained. 

8.1 The Model 

We use a core model defined by replacing the Nielsen-Olesen profile functions 

by, 

f{r) = u{r) = e{r-R,), (66) 

where 9{x) is the Heaviside step function. Clearly these profiles do not satisfy the 
Nielsen-Olesen equations, however it is still a perfectly valid exercise to investigate 
fermions in this background. The behaviour of the core model profiles for r > 
is identical to the asymptotic behaviour of the Nielsen-Olesen profiles, particles 
will therefore have their usual vacuum masses for r > Rc. While for r < Rc, both 
the Higgs and gauge fields vanish and all particles are massless. 

The scattering of fermions off the core model string has been investigated in 
and [^, for completeness some of their results will be rederived. 



8.2 Core Model Spinors 

To find the global solutions of the Dirac equation we solve the Dirac equations 
in the regions r < R^ (internal) and r > R^ (external) then impose continuity of 
the spinors at r = R^. For r > Rc the general continuum {E'^ > g^) solution is 
given by (^) with, z/ = n — A + a^. We now drop the subscript A in the definition 
of V for clarity. Similarly, the general external bound state solution {E'^ < g"^) is 
given by (|T]). 

The general regular internal solution reads. 



r<Rc 



a 



( UEr) \ 
-Jn+i{Er 










(67) 



For r < Rc the fermions are massless, the mass-shell condition reads E"^ = k"^ 
and the argument of these solutions is Er. These solutions are therefore suitable 
for use in the bound state energy region as well as the continuum regions. 

The spinor profiles must be continuous at r = Rc. For continuum states this 
gives matching conditions which allow us to determine the scattering coefficients 
analytically. 
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In the case of bound states we only have two regular asymptotic spinors and 
the system of constraints is over specified. Non-trivial solutions to the matching 
conditions are then possible if and only if, 



EK,+i{kR,)Jn{ER,) + 
kK,ikR,)Jn+i{ER,) 

gK,+iikR,)Jn-xiER,) 



gK,+i{kR,)Jn{ER,) 

EK,+^ikR,)J„.xiER,) + 
kK,{kR,).h^x+i{ER,) 



0. (68) 



This condition clearly puts restrictions on the energy values for which physical 
bound state wavefunctions exist and so generates a discrete bound state spectrum. 

Just as in the case of the actual Z-string there are also zero mode solutions. 
These occur in the = mode for the up quark and in the n = — 1 mode for the 
down quark. The internal solutions in the zero mode case are simply constants, 
while the external solutions decay exponentially at large r. 

8.3 Core Model Discretisation, State Pairing and Spectral Flow 

In appendix A we discuss in detail how to discretise the continuum by im- 
posing boundary conditions on the spinor profiles at some large radial distance, 
r = Ro- The advantage of the core model is that analytic expressions are known 
for the scattering coefficients and so an analytic expression can be written down 
for the discretisation condition. 

Each spinor in the discretised continuum consists of a linear combination of 
the two regular internal solutions. We denote the relative amount of each regular 
internal solution in the spinor by 7. In deriving the discretisation condition we 



obtain a quadratic equation for 7 ( |1471 ). By direct substitution of the scattering 
coefficients into ( |148| ), ( |149| ) and ( |150| ) we find explicit forms for the coefficients 
and the quadratic reduces to. 



'y'^igJn+lJn) + ^E{Jn+lJn-X " JnJn- 



gJn-X+lJn- 



0. 



We have solutions. 
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A+l) ± 



'^QJnJn 



+1 



where the discriminant. A, is given by. 



A — E'^{Jn+lJn-X — Jn-Ji 



n-A+lj 



X+l-Jn- 



(69) 



(70) 



(71) 



The argument of all the Bessel functions in (^)-(|7T|) is ERc. There are four 
cases to be considered, corresponding to | z/ |= ±1^ and to the two roots of the 
quadratic equation. Explicitly the discretisation condition becomes. 



J\u\{kRo 
N\^\{kRo 



for I V \ = C^v 



(72) 
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where, 

^± _ E {Jn+lJn-X + Jn-X+lJfi) J\u\ T J\v\\f^ — 2C,k J nJ n- \J .^gx 
~ E{Jr,+lJn-X + Jn-X+lJn)N\,\TN\,\^/E-2CkXJn-xN\,\+^ 

Throughout, unless the argument is exphcitly given, the Bessel functions of inte- 
gral order, associated with internal solutions, have argument ERf.., whereas those 
of orders involving | v \ have argument kRc- 

Using the small argument expansions of the Bessel functions it is fairly straight- 
forward to show that, 

X^\^\{k) tan(| i/ I tt) -1< i/ < 
XZ\^\{k)''^^Q -l<i/<0 
X%^{k)'7Z'Q v<-l. 

Similarly, using large argument expansions of the Bessel functions we find as 
k — > oo, 

K\ - cot J(n- I u I), Xr, - cot J(n- | u \ -A), (74) 



Xt\,\ ~ -cot|(|i/|-n), Xl|^|~-cot|(|i/| -n-A). (75) 



It is important to note that the infinite momentum behaviour of these functions 
is independent of the core radius Re- 

We can now calculate the fermionic energy difference between two core model 
strings with different core radii. This model calculation will illustrate many of 
the features of the full calculation, but in a more managable context. 

The eigenvalues are expected to be continuous functions of the core radius. Re- 
As a result a 1 : 1 correspondence should exist between the fermionic eigenvalues 
of the two strings. We first verify this by analysing a specific example. This 
provides insight into the physical intepretation of the singularity structure of the 
functions X^{k). 

As the core radius, i?c, is varied a continuous flow of the energy eigenvalues 
is expected, with the total number of fermionic states remaining constant. If a 
continuum state close to threshold loses energy and falls into the bound state 
region, a new bound state solution appears. Consequently, one of the infinitely 
many continuum states has become a bound state. It is essential when comparing 
two systems that the correct number of states are counted in each system so that 
the total number counted is the same in both cases, failure to do so would result 
in the occurrence of divergences which arc not removed by the counter term. 

For finite Rg the solutions of the discretisation condition are the points where 
the graph of ^^"'^(fc^"') intercepts the graphs of The function ^^"'^^fe^") behaves 
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like cot {kRo — — j) for kRo ^ 1, so there are an infinite number of branches 
with a width proportional to The singularities and zeros of clearly 
occur at the zeros of N\^\{kRo) and J\u\{kRo) respectively. 




Figure 1: Singularities in X^{k) and extra/missing states. The solid line shows 
a Type 1 and a Type 2 singularity in X^{k) respectively. States are found at the 
intesections of this line and the dashed line. 



The non-zero roots of J\y\{kRo) act as a suitable set with which the solutions 
of the discretisation condition ([7^ ) can be matched, and thus provide a means of 
counting the continuum states. Complications arise, as discussed in appendix A, 
due to the singularities in X^{k). For each Type 1 singularity there is one zero 
of J\y\{kRo) for which there is no solution of the discretisation condition [fig.|ll]. 
Conversely, for each Type 2 singularity there is one zero of J\^\{kRo) for which 
there are two solutions of the discretisation condition. The number of states will 
be denoted by (1 : l)|y| — q, where the notation signifies that each solution of the 
discretisation condition can be paired with a non-zero root of J\,,\{kRo) with q 
non-zero roots of J\y\{kRo) being left unpaired. 

This process is independent of the specific value of Ro- As Ro is increased the 
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density of branches of ^^^^'^^^^^"^ increases but the resulting state count does not 

change, because the functions X^^^^{k) are independent of Ro- 

To illustrate this we consider two core model strings, one with Rc = 0.5 and 
the other with R^ = 1. The parameters used are, 6^ = 0.5, g = 1.0. Summing 
over the ± solutions for modes n = —2. .2 for both core radii, we find 2(1 : l)|j,| 
states in all cases except for the function X~^{k) in the n = —1 mode. This 
difference arises because of the existence of a Type 1 singularity for Rc = 1.0 in 
X"'"(A;)[fig.0]. This leaves one non-zero root of the appropriate Bessel function 
J\^\{kRo) unpaired in the counting process. 



2 - 

Positive. 




Negative. 



12 3 4 5 

k 

Figure 2: X^{k) for the up quark with n = —1 and Rc = 1. 

To find the total number of states we must also include any massive bound 
states and zero modes. For Rc = 1.0 there is just one massive up quark bound 
state and it resides in the n = —1 mode. No massive up quark bound states exist 
for Rc = 0.5. Thus spectral flow has taken a continuum state and produced a 
bound state. 

The inclusion of the zero modes is more subtle as they only move in one 
direction along the string. A zero mode should only contribute | to the state 
count of the mode in which it arises. In the 2 + ID context the zero mode should 
therefore count as only half a state. This is an example of fermion number 
fractionization [^, ^ , but it has been motivated by looking at the full 3 + ID 
system. 

As expected, the total state count is the same within each mode hence the 
total number of string states is the same. It is worth noting that the state 
counting matches within each mode because there is no spectral flow between 
different angular momentum sectors as the core radius, Rc, is varied. 

The state counting procedure can be repeated using exactly the same tech- 
niques for the down quark. Again we find the same total state count in both 
backgrounds. 
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Having verified that the total number of fermionic states within a given mode 
does not change as Rc varies, we attempt to compute the fermionic energy dif- 
ference of the two core model strings. 

8.4 The Core Model Fermionic Energy Shift 

Using (|56|), the difference in fermionic energy within a given mode for two 
systems denoted 1 and 2 is, 

. ™. X ("Mm - "Vm) 1 , 1 dt _tE^ . ^ _is2 

AE (t) = ■ ■ / ^ > e — > e "•'2 

^ ' 871 T 4:71 Jr ^ 

dt 
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where the continuum contribution is encapsulated in the functions A{k) defined 
by, 



I cont 



A{k) = (3^71 + + cot-\X{k)) 
The 61, terms cancel and the continuum contribution becomes. 



1 r^e-^^Ekoo-/3r.)+-rci^^e--^f 
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2 r...„-*.^ f cot-\xm) 



I cont 



.(76) 



In the core model the large momentum limit of the functions X{k) is a con- 
stant which is independent of the core radius. Thus the cot^^(Xf''(A;)) 1^°"^* 
-cot-H^i^(^)) I™''* term tends to zero for large momentum. 

As discussed in appendix A, the coefficients of the first integral in (|77| ) can 
be deduced by simply counting the singularities of the functions X'^i^k), 

.± _ ± _ / No. Type 1 of X^^ X _ / No. Type 2 of Xt \ 
^''^ ^''^ " I - No. Type 1 of X^^ J U No. Type 2 of X^^ I 

For the up quark, the only singularity in the functions X^{k) occurs for the 
+ sign case in the n = — 1 mode. Thus we have, (3^^^ — ^2,00 = 0; cases 
except for the + sign case in the n = —1 mode, where the type 1 singularity in 
Xi+'"="\fc) implies that pt,^ - = 1. 

In all case except for the -|- sign case in the n = —1 mode, the energy difference 
for each solution within a given mode reads, 

A^^(r) = -— I -e-*^ dkke-'" , (78) 

where, 

F^'^ik) = (cot-i(X^"(A;) I,) - cot-i(X2±'"(A;) |p)) . (79) 
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In the special case [+ sign, n = 



— 1] we have, 



AE1=-\t) = — r% [e-*^--M - e-*^'[l - - r dkkF+'-\k)] 
Atx Jt t L tt jo 



2t 

IT Jo 



(80) 



The first term in (|80D is the standard bound state contribution due to the exis- 
tence of the massive bound state in this mode when Rc = 1.0. The second term is 
more subtle, it is due to the Type 1 singularity in X^'^^{k) which gives, for this 
particular case, f3i^^'^ — f32,oo'^ = 1 introduces an extra term from ([76|) . This 
term looks exactly like a massive bound state contribution except the state has 
threshold energy and the term arises with an overall minus sign relative to the 
massive bound state. In the state counting the Type 1 singularity in X+'~^(/c) 
ensured that there are the same number of states in each mode for both string 
backgrounds. Spectral flow considerations suggest that the natural state to pair 
up with a massive bound state is a continuum state with threshold energy. Such 
a state would make exactly this contribution to the fermionic energy. These ob- 
servations show that our formulation of the discretisation condition produces the 
desired state counting and matches states up in a consistent manner. 

Whilst an individual bound state contributes a manifest quadratic divergence, 
the coefficient of this divergence is independent of the energy of the state. Thus 
this pairing of a massive bound state with a threshold state cancels the quadratic 
divergences. This is to be expected from our naive arguments using a momen- 
tum cut-off: while the total fermionic energy must be quadratically divergent, 
the bound states only move along the direction of the string and so they only 
contribute a logarithmic divergence. 

The total continuum contribution made by a particular mode is determined 
by the sum F"*"'" and 

F'^(fc) = (81) 
± 

For all modes the function F"'(k) tends to zero at large k, while at small k we 
find -F"(0) = for all modes except for the n = — 1 mode, where F~^{0) = — tt 
(fig.0). This is an example of Levinson's Theorem which relates the number of 



bound states to the associated phase shift. In the notation used here, Levinson's 
Theorem states that the number of massive bound states of the first string minus 
the number of massive bound states of the second string is equal to — 

We should now subtract the counter terms to obtain a finite energy differ- 
ence. The quadratic counter term can then be compared to the behaviour of the 
fermionic energy as r tends towards zero and the quadratic renormalisation can 
be implemented. However there is a pathology with the logarithmic counter term 
which prevents a full renormalisation of the core model. We can however demon- 
strate several important features of the calculation, in particular the convergence 
with mode number. 
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Figure 3: F ^{k) for the up quark with 6yj = 0.5, g = 1.0 , i?c,i = 1-0 and 
Rc2 = 0.5. 



Using (p^), the divergent part of the quadratic counter term is, 

1 1 /""^ 1 
EQuadir) \m.= —- drrTi A[a^] = -^- drrUi{rf - f2{rf)^ (82) 
ioTT r JO 47r r jo 

where the subscripts refer to strings 1 and 2. Using the definition of the core 
model profile functions (|66D , we have. 



-Eqw(t) \Div= ^ — -Rci - -Rc2 = ^ — > (83) 
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for the example considered (i?ci = 1-0, Rc2 = 0.5). 

The coefficient of - in (|83D should be compared with. 



limr E Ai?"(r), (84) 



r^O 

n=— oo 



where A£'"'(r) is defined by ( ffSD and (|80|). The convergence of the sum over n 
is of critical importance both fundamentally and practically. Each A£'"(r) must 
be calculated numerically and hence in practice the sum must be computed up to 
some finite value of n, Umax, and at some finite value of r. We will demonstrate 
that for a given value of r, providing the value of rimax is sufficiently large, the 
higher terms in the sum over n are indeed small and convergence is ensured. 

The reason for this convergence can be seen by looking at the form of the 
energy shift functions, F^{k). Figure ^ shows the functions F'^{k) as a function 
of A; for n = 5, 10, 15 for the up quark, the features are however generic. As n 
increases, the functions F'^{k) remain close to zero for a larger range of k before 
the oscillations begin. This is due to the spinor profiles. They are Bessel functions 
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k 



Figure 4: The functions F"'{k) for the up quark with = 5, 10, 15. 



and so power law suppressed until their argument is of the same magnitude as 



their order 21 



The proper time regulator introduces the exponential suppression factor e~*'^^ 
into the energy integral. As the absolute value of the mode number increases, 
the region where the functions F^{k) are small (ie. k < n) becomes larger and 
there is increased suppression of the overall contribution. Thus overall we expect 
a suppression of the form, e~'^^ . 

Physically the modes with high angular momentum do not penetrate the core 
of the string and as a result are not affected much by changes within the string 
core. 

The total fermionic contribution made by a given angular momentum mode is 
found by computing the momentum integrals (|78|) and (|80|) and summing over the 
two cases which arise for each mode. The fermionic energy is found to decrease 
as the absolute value of the mode number n increases, but more slowly as the 
cut-off r becomes smaller. 

From a numerical point of view the momentum integral must also be truncated 
at some finite value. This value must be sufficiently large so that the exponential 
suppression has taken effect. As the value of the cut-off r is taken to zero, 
the values of both k and n at which the exponential suppression takes effect 
become larger, consequently we must extend the range of integration in k and 
include more modes inorder for the sum over n to converge. For the purposes of 
numerical calculation, if the momentum integrand is required to be suppressed 
by a factor tol from order unity at the truncation point, the minimum value of 
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momentum at which the integral can be cut-off is given by, 



llogtol 



The maximum mode number used should be of the same order as kmin to ensure 
that the contributions made by higher n modes are also negligible. 

It is now possible to compare the coefficients of the quadratic divergence 
and the quadratic counter term. Figure || shows a plot of the total regularised 
fermionic energy multiplied by the cut-off function of r. The 400 lowest 

angular momentum modes where included in the sum. The horizontal line is the 
coefficient of the quadratric counter term which has the value 0.02984 for the 
parameters used. 




Total energy multipled by the cut-off. 



Coefficient of tfie quadratic counter term. 



0.0005 0.001 0.0015 0.002 0.0025 

Cut-off : tau. 



Figure 5: The total regularised fermionic energy multipled by the cut-off. 



It can be seen from figure ^ that tE{t) fermion is tending towards the coefficient 
of the quadratic counter term, shown by the horizontal line. Thus we have 
control of the quadratic divergences. The curvature in tE(t) fermion becomes more 
significant as r becomes small due to the effects of the logarithmic divergence. 

Unfortunately we cannot implement logarithmic renormalisation in the core 
model. The problem arises because the discontinuous nature of the core model 
profile functions implies that the energy per unit length of the string is infinite. 
The Higgs field gradient term, for example, gives a contribution, 

/•OO /"OO 

/ drr{f{r)f= drr6\r - R^), (86) 
Jo Jo 

which is divergent. 
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Gradient terms of an identical form arise in the expression for the logarithmic 
counter term (^), giving rise to the wavefunction renormalisation of both the 
Higgs and gauge fields. As a result the logarithmic counter term is not defined in 
the core model background and so logarithmic renormalisation cannot be carried 
out. 

The core model has provided useful insights into the problem and tested 
critical aspects of our formalism. It has provided insight into the process of 
counting the fermionic states in the string background, in particular we have 
demonstrated that the total number of fermionic states is the same for two string 
configurations with different profile functions. The physical significance of the 
singularities in the functions X^{k) has been demonstrated: a Type 1 singularity 
corresponds to the occurrence of a massive bound state. The technology of fixing 
the discretisation condition uniquely has been shown to ensure that the correct 
number of states are included in the energy sum. In particular, we have seen 
that a massive bound state becomes naturally paired with a state at the bottom 
of the continuum, in accord with spectral flow considerations. An example of 
Levinson's theorem has been observed, with being the difference in the 

number of massive bound states in the two systems. Finally, the core model has 
demonstrated how the convergence of the sum over angular momentum modes 
arises. This convergence is critical for the success of the method, as numerical 
computations only ever allow a finite number of modes to be considered. 



9. The Perturbed String 

We now turn to the calculation of the fermionic energy difference between 
an actual Z-string and a perturbed Z-string. This will tell us what effect a 
perturbation in the fields that cause the Z-string to be non-topological has on 
the fermionic energy. The perturbations we shall consider involve the upper 
component of the Higgs field becoming non-zero. We set 0+ = g{r), where g{r) 
is some radial profile function which vanishes as r — oo. 



9.1 The Perturbed String Dirac Equation 

The upper component of the Higgs field couples the up and down quarks, 
hence the Dirac equation will have 8 components. Explicitly the (2 + 1)D time 
independent Dirac equation reads. 



E 



Ul 

dR 



u,R 



-9u(, 





V 9n 



-gu4>o 









( Ur\ 


-9d(t)+ 













-9d(p*o 




dR 





-gd4>o 




) 


\dL 1 



(87) 



where all the covariant derivatives are defined in the background of the unper- 
turbed Z-string and each of the spinors ur^ul, dji, di has two components. 
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Using the standard Z-string ansatz and 0+ 
decomposition, 



/ Mi? \ 

Ul 

dR 
\dL J 



E 



g{r) together with the mode 



we find the following set of 8 first order, ordinary differential equations for the 
radial spinor profile functions. 
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[dr + 



r 



+ a 



r 



r 
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'Ul - guf{r)u\ + gug{r)d\, 
U2 - guf{r)uf - gdg{r)d^, 



Eu^ = i-dr + ^ + af-^)u[ - guf{r)u^ - gMr)dl 



Edf 



-dr 



^""^^K ai^)d^ - g,f{r)d\ - g,g{r)u[, 



Edi 



(dr 



n 



i>{r) 



- + «i? — 



)di - gdf{r)d^ - gdg{r)u^, 



Ed[ = (dr + 



(n + 2) 



+ a\ 



z/(r) 



)d^-gdfir)d^ + gugir)u^ 



Edi 



-dr + 



(n 



1) 



+ a 



r 



)di - gdfir)d^ + gug{r)u 



4. 



(89) 

(90) 
(91) 
(92) 
(93) 
(94) 
(95) 
(96) 



This set of equations can easily be seen to decouple into the profile equations 
for the up and down quark in the background of the pure Z-string, (|30D to (|33D , 
if g{r) = 0. As g{r) vanishes at large radial distances, the system decouples 
asymptotically into 2 systems, one for the up quark and one for the down quark. 
Thus the asymptotic solutions are given by (^OD and (|4l|) for continuum and 
bound states respectively. 

It is important to note that the orthogonal subspace spanned by spinors with 
mode number n in (^) is also spanned by unperturbed up quark spinors with 
mode number (n+l) and the down quark spinors with mode number n as defined 
by the mode decomposition (^91) . 
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9.2 Perturbed String Discretisation 

Of the 8 solutions of the Dirac profile equations ( p9D to (0), 4 are regular 
at short distance. A spinor describing a physical state must consist of a linear 
combination of these 4 regular short distance spinors and will match onto a unique 
linear combination of the 8 regular large distance solutions. Denoting the regular 
short distance spinors by Si where i = 1,2, 3, 4, we have, 









[ 







(97) 



where a^^^ = [ai , a2 , , ) are the projection coerfacients found by matchmg 
to the large distance solutions. The matrices Mu^d are defined by the asymptotic 
Z-string continuum solutions, 



/ J\y\{kRo 



u,d 







N\,\{kRo) 







T^J\u\±i{kRo) TlJH±i{kRo) TjN\,\±i{kRo) TlN\,\±i{kR, 



\ 



J\u\{kRo) 



N\.\{kRo) 



\ ±lJ\,\^i{kRo) ±f J|H±i(A;i?o) ±fiV|,|±i(A;i?o) ±^N\,\^^{kRo) j 



where the appropriate subscript must be placed on both v and k. 



(98) 



Vu = n^a\, Ud = n + l + a^, E^ = kl. + gl 



By numerically solving the profile equations (p9| ) to (pq), each of the regular 
short distance spinors. Si, can evolved out to some sufficiently large distance, R, 
and matched to the asymptotic solutions. If we collect the up and down quark 
components into individual 4 component spinors, Vu{R) and Vd{R) respectively, 
we have, 

Vu = Sud, Vd = Sdd, (99) 
where a"^ = (ai, ^2, as, 04) and the scattering data is encoded in the matrices, 

Su,dj 



Su,i 



I V'l,! 

/ u,d 
^2,1 

/ u,d 
^3,1 

/ u,d 



I u,d 
^1,2 

I u,d 
^'2,2 

/ u,d 
^3,2 

/ u,d 



I u,d 
V^l,3 

^2,3 

; U,d 

V^3,3 

; U,d 

^4,3 



; U,d 

^1,4 
/ u,d 
^2,4 

; U,d 

^3,4 

; U,d 

^44 



Here ip^f is the value of the i-th component of the up/down quark spinor at 
r = R found by evolving the j-th regular short distance solution. From (|97|) the 
4 component spinors Vu and Vd can be written as. 




Mu 
Md 




(100) 
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By eliminating a using (0), we can express a„ or Ed in terms of the other, 

ad = Mdu, (101) 

where, 

M = M^^SdS-'Mu. (102) 

Thus only 4 out of the 8 asymptotic projection coefficients are linearly indepen- 
dent. We arbitrarily treat the components of a„ as our arbitrary coefficients. 
These degrees of freedom are of course related to the components of a. In fact, 
there are only 3 degrees of freedom because the spinor is only defined up to an 
overall normalisation constant. 

As with the pure Z-string, the next step is to discretise the continuum by 
imposing boundary conditions on the spinor profiles at r = i?o. We double up 
the conditions used for the pure Z-string, this allows much of the technology 
developed in section 6 to be applied to the perturbed string. Specifically we set 
the first, third, fifth and seventh components of the 8 component spinor to zero 
at r = Ro. Clearly as the perturbation is switched off these conditions become 
identical to those used in the case of the unperturbed string. 

Applying the boundary conditions at Rq we find. 



J\y^\{KRo) _ al _ al 



and, 

J\vAkdRo) ai af 



(103) 



(104) 



Ni,,iikdRo) af af 

Individually ( |103| ) and (|104|) are identical in nature to the discretisation con- 
ditions which arose in the case of the pure Z-string. However ( p.03|) and ( |104| ) 
must be solved simultaneously due to the non-trivial coupling between the up 
and down quarks induced by the upper component of the Higgs field. 
We seek a discretisation condition of the form, 



J\u\{kuRo) 



-X{k), 



in order to employ the methology developed in section 6. From ( |101| ) we have. 



(105) 



af Mi,i^ + Mi,2 + Mi,3 Jf| + Mi,4| 
and 

nd M4.1 ^ + M4.2 + M4.,^^ + M^A 

(106) 



^4 



ai M2,i^ + M2,2 + M2,3^^ + M2,4^ ' 
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where has also been ehminated using ( |103| ). 

Implementing part of (|104|) by equating ( |105| ) with (|106|) , we find a quadratic 

Q^U qU 

equation for thus we can consider ^ to be a function of Substituting 
the solutions of this quadratic into ( |105| ) and ( p,06| ), the energy of the continuum 

qU 

states and the ratio ^ are the only remaining independent degrees of freedom. 
We must now impose the 2 remaining constraints, 

to determine the allowed energy E and the ratio 

A solution to this system has only been found for degenerate masses, Qu = Qd 
and hence fc„ = kd- For degenerate masses and kRo ^ 1, using the standard 
large distance expansions of the Bessel functions, we have. 



«3 

02 af 



coi{kRo - 0„) = --i, cot{kRo - e, 
where, 

Ou,d = I +1)- (108) 

Using standard trigonometric identities to eliminate kRo gives, 

;;^ + tan6'„ ;;| + tan6'd 

0. (109) 



^tan^„-l ^tan^d-1 



The only independent variable in ( |109| ) is the ratio ^ and we have a one dimen- 



sional root finding problem. This process, although conceptually straightforward, 
requires a significant amount of computational effort due to the singularities found 
in ( p.09| ). Numerically solving ( |109| ) for a range of continuum energy values pro- 
duces functions %{E) which can then be used in ( |103| ) to provide a discretisation 
condition of a similar nature to that of the pure Z-string. 

For a given value of the energy, 4 solutions of (|109|) are expected. Denoting 
these solutions by X^'^k) where a = 1,2,3,4, for kRo ^ 1 the discretisation 
condition reads, 

cot{kRo - e^) = -X"(A;), (110) 
which can be inverted using (^) to give, 

kiRo = 171 + p^n + 9u + cot-i(-X"(A;)) 1^*, (HI) 

with the value of given by ( |156| ). 
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Using the analysis of section 6, the continuum contribution to the fermionic 
energy of the perturbed string made by a given angular momentum mode is, 



Ro dt 



°o dt 



1 ' 
Air 



1 dt 



r\t o 



l-K Jr t2 

6„ dt 



je-'^ dkke-"" cot-^(-X"'"(A;)) |;;°"* . (112) 



As in the case of the unperturbed string, the first two terms cancel when two 
systems are compared. 

We must also consider the bound states of the perturbed string by looking 
for suitable linear combinations of the regular short distance spinors that match, 
at large radial distances, to spinor profiles that decay exponentially. In the per- 
turbed string there are 4 regular short distance spinors in each mode, we denote 
these by Si. The perturbations we are considering are localised to the string core 
so, just as in the case of the continuum, the 8 linearly independent large distance 
solutions are the same as those which arose in the case of the unperturbed string 

&■ 

As in the continuum, we have a doubling up of states and matching conditions. 
The existence of a massive bound state solution requires that the coefficients of 
the exponentially divergent solutions vanish when a suitable short distance spinor 
is evolved to large distance. Taking a general linear combination of the 4 regular 
short distance spinors gives the bound state condition. 



0, 



(113) 



where the summation convention is used and the asymptotic states are defined 



in (p|). These conditions lead to the determinant equation. 



"1,1 
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^1,2 

^1,1 



"1.2 



''2,1 



'-2,2 

^2,1 
d 



''3,1 
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''3,2 
''3,1 



''4,1 
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^,2 
^d 
^,1 



''3,2 "4,2 



0. 



(114) 



We find two bound states in the perturbed string background. 



E' 



n=-l 



0.0505, E' 



n=0 



0.9992. 



(115) 



n = —1 corresponds to the mode in which the up and down quark zero modes 
reside if the perturbation is switched off, thus, as first shown in [|r^, the per- 
turbation lifts the energy of the zero modes to form a low energy bound state. 
The index theorem 121111 which ensures that there are zero mode solutions in the 
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background of the pure Z-string no longer applies because the perturbation in 
the upper component of the Higgs field is exactly one of the fields which makes 
the electroweak string non-topological. In terms of the spectral flow of fermion 
energy levels, as the perturbation in the upper component of the Higgs field is 
turned on, the up and down quark zero modes, which are restricted to move in 
opposite directions along the length of the string, are combined to from a low 
energy massive bound state. The number of fermionic degrees of freedom asso- 
ciated with the zero modes is therefore unchanged. Similar behaviour has been 
observed in the context of GUT string zero modes mixing with massless modes 



Denoting the set of massive bound state solutions in a given mode by {En^i}, 
the total contribution made by the mode to the fermionic energy becomes, 

1 dt 



27r2 



/ -e"*^ / rfA;A;e-*'= Vcot-^(-X"'"(A;)) 1;^°"* . (116) 
Jt t Jo 



9.3 The Perturbed String Counter Term 

The calculation of the counter term for the perturbed string proceeds along 
the same lines as the corresponding calculation for the pure Z-string. The 4 
component Z-string Dirac Hamiltonian is replaced by the 8 component Dirac 



Hamiltonian defined in (0). The counter term is the sum of the individual 
up and down quark counter terms in the background of the pure Z-string (as 
defined by (Q), ( |65D and (|63D) together with some modification arising from the 
perturbation. 



By squaring the Hamiltonian in ( P7| ) and writing it in the form, 

= -{d, + r,)(9, + r,) + a(x), (117) 

we can apply the formalism for the heat kernel expansion developed in section 7. 
Explicitly we have, 

Ti = iqZi diag (0^/2x2, -alhx2, "^^2x2, -aihx-i) , (118) 

and. 
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where, 



au{x) 



-iqa\{diZi - a^a'diZj) 

-igu(T'di{f{r)e-'^) 
-g^qa^ZJ{r)e-''' 



) \ 



iguCT%{f{r)e'' 
-g^qa^ZJ{r)e^' 



iqal{diZi - a^a'diZj) 



ad{x) 



I -iqajiidiZi - a^a'diZj] 
Wj{r)^ + 9l9{r)^ 

-igd(T'di{f{r)e'^) 
+gdqa'ZJ{r)e'' 



and, 



ai{x) 



-i9d9{r)a'dig{r) 



tgdCT%{f{r)e-^') \ 
+gdq(y^Zif{r)e~''^ 

iqaiidiZi - a^a'diZj) 
+9lf{rf + 9l9{rf J 



9uq9{r)a'Zi{a\ + ai) \ 
-igucr'dig{r) 



\ -q9d9{r)a'Z,{aji + al) {gj - gl)gir)f{r)e~'^ J 



(121) 



(122) 



After a similar, but more lengthy, calculation to that presented in appendix B, 
we find the heat kernel coefficients in the perturbed string background are. 



and, 
Tr[a2] 



nai] = -<9l + 9l){f{rr + 9{rn 



+ {9if{rr+9ig{rrr+{gif{rr+gig{rrr 

+ {9lf{rf + 9l9{rff + {9lf{rf + 9l9{rff 

7(r)Kr) + l]V 



(123) 



+ '^{9i + 9i)\U{ryY + 



+ 2i9i-giygiryfir 



+ ^9i\i9ir)r + {al + < 



dx2 f ^ir)g{ry 



+ 2gi igiryy + ia% + al 



u\2 



'u{r)g{r) 



(124) 



We see that ( |123| ) and ( [124| ) reduce to the sum of the individual contributions 
made by the up and down quark in the background of the unperturbed string 
when the perturbation is set to zero. 
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Having found expressions for both the fermionic energy of the perturbed string 
and the corresponding counter term, these results can be combined with the 
corresponding calculations for the unperturbed string to compare the fermionic 
energy of a perturbed string to that of an unperturbed string. 

9.4 The Change in the Ground State Energy 

In order to obtain a well behaved mode decomposition of the energy difference, 
we must examine the angular dependence of the spinors. The space spanned by 
the nth modes of the perturbed string energy eigenstates (^) is also spanned by 
the set formed by the union of the (n + l)th up quark modes and the nth down 
quark modes in the background of the unperturbed string, as defined by (pQl) . 
Physically this means that the spectral flow caused by the perturbation only 
mixes the (n + l)th up quark mode and the nth down quark mode. Thus the 
energy difference is found by subtracting the sum of the contributions made by 
the (n+ l)th up quark mode and the nth down quark mode in the background of 
the unperturbed string from the contribution made by the nth perturbed string 
mode. Matching modes in this manner ensures that the fermionic states of the 
two systems can be put into a 1 : 1 correspondence. 

This coupling of modes is consistent with the mixing of the up quark zero 
mode (n = —1) and the down quark zero mode {n = 0). 

Given this coupling of the modes, it is convenient to make the following map- 
ping on the unperturbed up quark mode number which was used in section 4, 

n„ ^ n — 1. 

Again we split the energy difference into a bound state contribution and a con- 
tinuum contribution. For the bound states we have, 

^Kow(-) = + ^E(«.)^r(-i,r«j^) 

-i(|:/r.-E/^^5-E/3-)/r(-i,A), (125) 

where {E^j^}, {E!^ j^} and {E^j^} are the set of positive energy massive bound state 
solutions for the perturbed string, down quark and up quark respectively with 
mode number n . a^Af u ^^zm d denote the number of up and down quark zero 
mode solutions respectively in the unperturbed background with mode number n. 
The values of the integers /3oo required to fix the discretisation condition uniquely 
at large momentum, are determined by counting the singularities in the functions 
X{k) and using ( |156| ). 
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The continuum contribution to the energy difference is given by, 



27r' 



^ . r^e-*^' rdkke-''" (^cot-i(-X"'"(A;)) |^°"* 

Jt t Jo \a=l 



-^cot-i(-X"''^(A:)) 1;;°"* -^cot-^(-X"'"(A;)) |;;""*) .(126) 
± ± J 

The second term in (11261) is of a famihar form, while the first is a new feature of the 



perturbed string. Treating ^ as the independent variable gives a discretisation 

°'2 

condition of the form in ( [L 1 1| ) . This has the same form as the discretisation 
condition for the up quark in the unperturbed string background. Had we used 
as the independent variable instead, the discretisation condition would have 
taken the form, 



kiRo = ZTT + /^ooTT + + COt-^(- X (k)) 



cont 

\p 5 



which is of the same form as the discretisation condition of the down quark in 
the unperturbed string background. Thus, by working with we are taking 
the up quark states of the unperturbed string as our reference points in both 
sectors. The first term in ( |126| ) represents the shifts of the down quark states 
of the unperturbed string relative to these up quark reference points. There 
are of course corresponding factors within the perturbed string phase shifts, so 
that overall the energy difference is independent of the reference points used. It is 
most natural to absorb this term into the integral over the continuum momentum. 
Using ( |10^ ) it is easily seen that, 

0..-0., = ^i\u.,\-\u,\), (127) 

and so, 

1 i-OO 

5^cont = I dkkm T{e + g'))F-ik), (128) 

where, 

F'-ik) = Ecot-i(-X^"„(fc))|-*-^cot-^(-Xj„(A:))|-"* 

-T.^ot-\-X^Jk)) |-"* -vrd I - I I). (129) 

± 

After some simplification, the counter term for the energy difference is given 



by, 



EcT{T,fi) = Hdrr ({- - -)AK] +log(^)A[a2]') , (130) 
167r Jo \ T jJ, T J 
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where, 

^[a,] = -A{gl + gl)g{r)\ (131) 

and, 

AH = 2{gt + d%g{rf + 2f{rfg{rf) + 2{gl + gl){g\r)f 



+2 {glial + air + aliai + al?) (^^^M^ j . (132) 

The degenerate mass case under consideration is obtained by simply putting 
g^ = g, = g in (^) and (Q. 

The counter term in ( [1301 ) still contains the arbitrary scale /i. As discussed 
in section 7, the total energy is independent of the particular value of /z because 
the parameters in the theory also depend upon /i in exactly the manner to ensure 
that physical quantities are independent of fi. The value of /i is set by the typical 
energy scale in the problem, which we take to be the mass of the Z boson. By 
fixing the vacuum expectation value of the Higgs field to he u = 247 GeV, we 
have, 

76.9 GeV , , 

I sm26'^ I 

Taking 6^ = 0.5, we have ~ 0.311 in units of the vacuum expectation value. 



2 



giving. 



/X = ^ ^ 10.32. (134) 



Finally, performing the integral in (|130|) for the profile functions shown in figure 
1^, we find that the counter term is given by, 

EcAr) = ai f r - 10.34) + a2 log , (135) 



T 



where, 

ai = -3.97886 x 10~^ 02 = 5.76546 x 10"^ (136) 
The specific parameter values used in the calculation were, 

6^ = 0.5, g, = g^ = 1.0, p = ^ = 0M3. (137) 

The actual background profile are shown in figure |^. The Higgs and gauge field 
profile equations /(r) and z/(r) were found by solving the Nielsen Olesen equation 



(|I0|), (|TT]) for the parameters given above using numerical relaxation The 
perturbation in the upper component of the Higgs field was taken to be a standard 
gaussian function, so that it is non-zero inside the core of the string but decays 
rapidly outside the core; 

g{r) = 0.1e-'\ (138) 
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1 2 3 4 5 

Figure 6: The background profile functions used. 



As discussed above, we find two massive bound state solutions in the back- 
ground on the perturbed string. The energies of these bound states are found to 
be, 

E^^Q = 0.99918, E^^_^ = 0.0505. (139) 

In the background of the unperturbed string the only massive bound state is a 
down quark bound state with n = and energy, 

E^^Q = 0.99918. (140) 

This energy is, to this accuracy, identical to that of the massive bound state in 
the n = mode in the perturbed string background. Hence these two state are 
naturally paired in the bound state energy contribution. 

The zero modes that appear in the pure Z-string background naturally pair 
up with the low energy, n = —1, massive bound state in the perturbed string 
background. 

Given this natural pairing of the bound states, all the other contributions in 
( |125|) should vanish for each n. That is. 



E /^^p - E /^S - E = 0, (141) 



jn,± 
oo,d 

a=l ± ± 



where the values of the integers (3 are determined by the singularities of the 
functions X"-{k) through (|156|) . This is indeed found to be the case. Physically 
this means that in the spectral flow between the unperturbed string and the 
perturbed string, no continuum states drop below threshold energy to become 
massive bound states. In accordance with Levinson's theorem, F"'{k = 0) = for 
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all values of n, as there is no difference in the number of bound states between 
the two systems. 

The functions F^{k) that define the shift in the energy of the continuum states 
are defined by (|129|) . They have been found numerically using the formalism 
developed in sections 3-6. Figure shows plots of the functions F"(/c) for n = 
0, 1, 2, 3, 4, 5, similar behaviour is found for negative n. 




There are a number of features of the functions F"'{k) worth noting: Except 
for a small range around zero momentum, the functions F^{k) are positive for 
all values of n. All other factors comprising the integrand for the continuum 
contribution for each mode are positive definite. However there is an overall 
minus sign, so each mode contributes a negative definite amount to the regularised 
fermionic energy. 

The values of = 0) are indeed zero for all n (closer inspection of the 

behaviour of F^{k) for small k reveals the function does indeed vanish). 

The functions F^{k) are enveloped by the functions with smaller absolute 
mode number. As all the other factors in the integrand for the energy contribution 
are independent of n, this implies that as the absolute value of the mode number 
increases the contribution made by the modes to the regularised continuum energy 
decrease. 

In units of the Higgs VEV, F^{k) ~ for k less than about the absolute value 
of n. This behaviour was noted in the case of the core model and is the key to 
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the convergence with n. Physically, as the spinors with high angular momentum 
are small within the core of the string for this range of momenta, they do not 
couple strongly to the perturbation. 

The functions die away as k becomes large. Numerical fitting to a 

function of the form ^ + fl + ■ ■ ■ suggests a | fall off. 

We can now perform the integrals in ( |128|) for specific values of the cut-off 
r and obtain the regularised fermionic energy. Numerically, the most intensive 
part of the calculation is the production of the functions F^{k). These functions 
were tabulated for n between —20 and 20. 

Figure ^ shows the continuum contribution to the fermionic energy made by 
each angular mode for n between —10 and 10 and for three values of the regulator 
r. We see that the contribution made by each mode rapidly decreases as the 
absolute value of the mode number increases, this implies rapid convergence of 
the the sum over n. However, as r — 0, the contributions from the low n modes 
can be seen to get larger, the fall off as \n\ gets larger is slower, and more terms 
must be included in the sum over n, as discussed in section 8.4. 
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Figure 8: A plot of SE^^^^ir) as a function of n for r = 0.2, 0.1, 0.05. 




We find the regularised fermionic energy for a given value of r by adding 
together the contributions made by the bound states and the continuum in each 
mode and then summing over a sufficient number of mode to ensure convergence. 
The renormalised fermionic energy is then obtained by subtracting the counter 
term. 

There are several sources of numerical error in this calculation. These arise 
due to truncation of the sum over angular modes, truncation of the momentum 
integral, finite tabulation density of the functions and numerical errors in 

the values of these functions. 

The first two errors limit how small the cut off r can be taken. Figure ^ shows 
the regularised continuum energy scaled by the cut off as a function of the cut 
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off for three different values of the maximum absolute value of the mode number 
taken in the sum. As the number of modes included in the sum decreases, the 
continuum contribution to the regularised fermionic energy becomes less negative. 
As the contribution made by each mode to the regularised energy is negative 
definite, the truncation in the sum over n induces a systematic error in the total 
regularised fermionic energy. As the number of modes included increases, r can 
be reduced further before significant errors arise. 




Figure 9: Plots of TdEamti'T') as a function of r for different truncations in the 
sum over the mode number n . 

As discussed in section 8.4, the combination of the regulating factor, e~**^^, 
and the fact that the functions F^\k) arc small for k < n, leads to an exponential 
suppression of the contribution from high angular modes by a factor of the form, 
g-rn Pqj. ^ maximum mode number of 20 and a suppression factor of order tol, 
the smallest value of r which can be used with confidence is given by, 

r ^ (142) 
-400 ^ ^ 

To quantify the errors arising from the intcrgal over the momentum wc varied 
the interpolation procedure and the tabulation density. Second and third order 
polynomials were used to interpolate between the tabulated points and a second 
order polynomial fit was made to only half the data points. The error was then 
estimated from the maximum absolute difference between these three methods. 
Again the error increases as the cut-off is taken away. 

The tabulation density used, the truncation k„iax = 50 and the truncation 
nmax = 20 all lead to errors of order 10~^ at r ~ 0.04. For smaller r values the 
errors grow rapidly. 

Using a quadratic fit to the nodes with r > 0.04 and extrapolating to r = 0, 
we find, 

SEfermion = -0.0033 ± 0.00002, (143) 
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where the error has been estimated by fitting to fewer nodes. 

We have an exphcit calculation of the difference between the fermionic ground 
state energies of the perturbed string and the unperturbed string. This energy 
difference is negative, as conjectured by Nacuhch |jl2|. The filled Dirac sea further 



destabilises the string with respect to a perturbation of this type. 

9.5 Scaling With Perturbation Size 

Repeating this calculation for a perturbation that is five times larger than the 
previous one, we find the same number of bound states and that the functions 
F^{k) have the same characteristic features but larger typical values. The typical 
increase in these values is of the same order as the increase in the perturbation 
squared (i.e. ~ 25). 

Figure [1^ shows a plot of the renormalised energy at a function of the cut-off. 
The extrapolation to obtain the value of the fermionic energy was again carried 
using a quadratic fit. The renormalised fermionic energy is found to be, 

5Ef ermion = -0.0667 ± 0.0002. (144) 



.U -0.02 - 



Figure 10: A plot of the renormalised fermionic energy 6Eren{'^) for the larger 
perturbation. 



Again the fermionic energy can be seen to be negative and hence does not 
stabilize the Z-string. 

Comparing the magnitudes of the two fermionic energy shifts, we see that the 
effect has increased by roughly the same order as the increase in the perturbation 
squared (i.e. ~ 25). This is consistent with the suggestions of Naculich [Q. 

We have calculated the renormalised fermionic energy difference between the 
perturbed and pure Z-string. In both cases the fermionic energy is found to be 
negative and hence has a destabilizing effect on the string. Our values appear to 
scale roughly with the square of the magnitude of the perturbation. 
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9.6 Populating Positive Energy States 

Having calculated the energy shift of the fermionic vacuum, we can now con- 
sider the effects of populating positive energy bound states. The inclusion of 
positive energy bound states is fairly straightforward, the positive energy states 
are filled up to some given value of the momentum along the length of the string, 
k. As the perturbation increases, the energy of these states tends to increase and 
the extra contribution to the fermionic energy shift will be positive. At some 
value of k the contributions from the fermionic ground state and the populated 
positive energy bound states cancel. In this situation the fermionic energy has 
a neutral effect on the stability of the string. For higher values of k the overall 
energy shift is positive and the fermions enhance the stability of the string. 

It was shown in section 3 that the fermionic contribution to the energy of a 
field configuration is given by, 

J2i-Er + nrEr) , 
r 

where the sum is over all the positive energy states and rir is the occupation 
number of the positive energy states. If an adiabatic perturbation were turned 
on then the change in the fermionic energy is given by, 

where (5-E<jea is the sum of the shifts in the negative energy states computed above, 
and SEpos is the shift in the energy of the occupied positive energy states. For 
an adiabatic change, the occupation numbers of the states remain constant and 
we have, 

SEpos = '^nr {E'j. — Er) , 

r 

where E'^. is the shifted value of Er when the perturbation is switched on. In the 
spectral flow picture these state flow into one another. 

In the case of the string the states of interest are the zero modes in the 
unperturbed string background and the low energy massive bound state which 
the up and down quark zero modes combine to become when the perturbation is 
turned on. These are the only states between which there is a significant energy 
shift. 

We start in the pure Z-string background, with up quark zero mode states 
with momentum in the range [0, k] filled, and down quark zero mode states with 
momentum in the range [—k, 0] filled. In the perturbed string this gives low, 
positive energy bound states with an absolute momentum up to k filled. Figure 
Tl| illustrates the population of positive states for the perturbed string (low energy 



massive bound state) and the unperturbed string (up and down quark zero modes) 
up to a value of k = 2.5. If the low energy bound state has mass 6m, then 
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Figure 11: A typical scenario illustrating the population of positive energy states. 



a state with momentum along the string satisfies the mass-shell condition, 
E"^ = k1 + 5m? . The zero modes move at the speed of light along the length of 
the string and therefore their energy, satisfies, E = kz, as illustrated by the 



45° zero mode lines in figure 11 



Discretising the momentum along the string as in section 5, we have. 



5E, 



pos 



E 



k?, + 5m2— I kr, 



(145) 



where, kr, 



2Tm 
L ' 



and L is the length of the cylindrical box we are considering. 



Taking L to infinity and defining the maximum momentum k by, k 
we obtain, 

1 rk 



6E, 



pos 



TT JO 



dk{VWT5^- \k\) 



6m 



2 r 



1 + 



k^ 



Sm'^ 



+ 



k ^ 
5m 



1 + 



A;2 



_k_. 

6m' 



6Epos is positive as expected, thus populated positive energy states help to 
stabilise the string. 

There is also a shift in the energy of any massive bound states when the string 
is perturbed. The energy shift for these states is simply found by considering 
the (possibly fictitious) zero mode as a reference and taking the difference of 
calculated with the initial and final masses. 

It may be possible to introduce a sufficient population of positive energy 
bound states and zero modes to stabilize the string against the negative curvature 
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induced by the Dirac sea and the bosonic instabihties. This question is considered 
in ref-H- 



10. Conclusions 

We have systematically calculated the order h fermionic energy shift when an 
electroweak string is perturbed. This calculation is the first to take into account 
the effects of the Dirac Sea. Our approach allows the infinite volume limit to 
be recovered analytically. We have developed a state counting procedure that 
correctly counts states flowing from the continuum to bound modes and allows 
the fermionic energy to be renormalised. The momentum shift functions we 
evaluate satisfy Levinson's theorem. Convergence of our numerical sums and the 
finiteness of the renormalised energy shifts have both been demonstrated. 

The Dirac Sea energy shift when a pure Z-string is perturbed is found to be 
negative and scales roughly as the size of the perturbation squared. The Dirac 
Sea thus provides a term with negative curvature to the effective potential and 
has a destabilising effect on the string. 

However, populated positive energy states have their energy increased by the 
perturbation and hence have a stabilising effect on the string. The total fermionic 
effect may vanish if there is a sufficient population of positive energy bound states 
on the string. A larger population of bound states leads to the fermions having 
a stabilising effect overall. This raises the facinating possibility that electroweak 
strings may be stable if they carry a fermionic current. 
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Appendix A: Discretisation Of The Continuum 

To carry out the sum over the continuum states we first impose conditions 
on the spinor profile functions on the surface of a cylinder of radius Rq centred 
on the core of the string. The conditions used select a discrete subset of the 
continuum states and we trace over this discrete set of states. The full continuum 
is recovered as Ro ^ oo. This can be done analytically, so that the actual trace 
over the continuum can then be done in the physical, infinite volume limit. The 
use of the discretisation condition is merely an intermediate step in setting up 
the sum over the continuum states, it is critical as it ensures that the correct 
density of states is used. 

To discretise the continuum we demand that the 1st and 3rd components of 
the spinors vanish at some large radial distance r = Ro. The physical answer 
clearly should not depend in anyway upon the discretisation condition used and 
alternatives could have been used. The critical point is that as Ro is taken to 
infinity, the full continuum should be recovered, so all states are included in the 
sum. That this is true for the condition used here will be become clear shortly. 

The Dirac equation in the Z-string background has two regular short distant 
solutions and two irregular solutions. Denoting the regular short distance spinors 
by Si where i = 1,2 then each of these will asymptotically match onto a unique 
linear combination of the asymptotic solutions given in (^OD, 



S,. 



( J|.|(/t) N\,\{k) \ / a,,i \ 

T^J\u\±l{^^) T^J\v\±i{t^) TfA^|H±i('«) TfA^|H±i('«) «*,2 

J\v\{i^) Ny{K) ai^3 

V ±1 J|,|±i(/s:) ±f J|,|±i(k) ±1 Ar|,|±i(/s:) ±fN^,^±^iK) J \ a,,4 / 

where, k = KRq, the upper sign corresponds to the case \ u \= u and the lower 
sign I z/ 1= — z/. 

By taking an arbitrary linear combination of 5*1 and 5*2 and setting the 1st 
and 3rd components of the spinor to zero at r = Ro, we obtain the constraints, 

a{J\u\{kRo)ai^i + N\^\{kRo)ai^3) + P{J\u\{kRo)a2,i + N\y\{kRo)a2;3) = 0, 
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a{J\,y\{kRo)ai^2 + Ni^\{kRo)ai^4) + (3{J\^\{kRo)a2,2 + N\y\{kRo)a2,i) = 0, 

from which it is easily seen that, 

J\y\{kRo) _ (aoi^a + /5a2,3) _ (aai_4 + /?a2,4) 
A''|,,|(A;_Ro) (aai,i + /5a2,i) (aai,2 + /?a2,2) ' 



(146) 



Only the ratio 7 = ^ is of interest and it is straightforward to derive the 
following quadratic equation for 7, 

A72 + E7 + C = 0, (147) 

where, 

A = ai,3ai^2 - cti,icn,4, (148) 

B = 01,302,2 + 02,301,2 - 01,102,4 - 02,iOi,4, (149) 

C = 02,3^2,2 - 0-2,102,4- (150) 

The solution to ( |14?| ) can be substituted into ( [L46| ) to produce two solutions for 
each particle, 

-X^{k), (151) 



J\u\{kRo) 
N\,\{kRo) 



where. 



{-B ± - 4AC)oi,i + 2Ao2,i 

The left hand side of the discretisation condition is approximately cot{kRo), 

thus for finite Ro there a discrete set of string states {fc,} with i G {1, 2, 3, } 

and spacing proportional to The density of the solutions becomes infinite as 
Ro ^ 00 and the full continuum is recovered. 

It is possible to calculate the solutions of (0) numerically, for a fixed value of 
Ro and then compute the fermionic trace. The process can then be repeated for 
different values of Ro and the limit Ro —>■ 00 taken numerically This is not 
the approach used here, instead the limit will be taken analytically. This saves 
the computation being repeated for different values of Ro. 

For momentum ki such that fcj-Ro 3> 1, we can take the leading order asymp- 
totic behaviour of the Bessel functions and write the discretisation condition as, 

cot{kiRo-e,) = -X{h), (153) 

where. 

This can be inverted to give, 

kiRo = m{i)7i + e^ + cot-\~X{ki)) \p, (154) 
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where \p denotes the principle range of the inverse cotangent, throughout this 
is taken to be [0,7r]. m{i) is some integer valued function which specifies which 
branch of the cot function the solution lies on. 

It is found that X{k) is singular for certain values of k, the physical signifi- 
cance of such singularities is discussed in sections 8 and 9. A singularity in X[k) 
corresponds to a discontinuity of tt in the function cot^^(— |p. As we take 
the limit the Ro oo, the discrete set of solutions, {fcj}, tends towards a contin- 
uous spectrum. This continuum spectrum, as its name implies, is continuous in 
the momentum, thus the difference between consecutive elements of the set {fcj} 
should tend to zero as the limit is taken. The requirement of continuity of the 
spectrum and the existence of singularities in X{k) imply that some extra factor 
must be included. This extra function, m{i), must jump in integer steps around 
the singularities of X{k) to ensure that the spectrum defined by {ki} is indeed 
continuous as Ro is taken to oo. 

It is convenient to write, 

m{i) = i + I3{i), 

as this allows an evenly spaced free momentum k" to be defined, 

k^Ro = in. 

The function f3{i) will be combined with the cot^^(— X^(A;j)) \p term to define a 
continuous function. 

In regions where X{k) is finite, we find one solution per branch of the cot 
function, thus l3{i) is a constant. However, if X{k) is singular, there can be 
branches of the cot function with no solutions or two solutions (fig.|l]). We denote 
these singularities as Type 1 and Type 2 respectively. 

Each branch of the function jg bounded by two consecutive roots of 

N\y\{kRo). Each solution of the discretisation condition lies along one of these 
branches and is thus bounded by two consecutive roots of N\j,\{kRo). Denoting 
the ith root of N^{kRo) by aj', the ith solution of ( pi]) will satisfy the inequality, 

<+5 < kiRo < (155) 

where 5 is an integer determined by the singularity structure of X{k). 

If X^{k) has no singularities, it is easy to see that the ith string state is 
bounded by the ith and the {i + l)th roots of N,j{kRo) giving 5 = 0. A Type 1 
singularity procludes a solution along one branch, thus above the singularity the 
ith solution is bounded by the (i -|- l)th and {i + 2)th roots of Ny{kRo)- In this 
case 5 = 1. Conversely, a Type 2 singularity gives two solutions on one branch. 
Thus above the singularity the ith solution is bounded by the {i — l)th and the 
ith roots of Ny{kRo), giving 6 = —1. 
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By counting the singularities of each type, we can determine the value of 6 at 
momenta greater than that of the highest momentum singularity of X{k), 



6 = ( No. Type 1 Singularities — No. Type 2 Singularities). 

The asymptotic value of /?, Poo, can then be set such that the solutions of the 
discretisation condition fall in the correct range. To leading order, the expansion 
for the large zeros of a Bessel function gives. 



rmr + 9^ — 71 z/ > 

mvr + — TT + int[—v\ < 



To place the high momentum solutions on the correct branches of the cot function 
we therefore set, 

/9oo = 5 - 1. (156) 

Having fixed the asymptotic value of /3, we can now consider the effect of 
singularities in X{k). Physically the occurrence of such singularities changes the 
number of continuum states which exist. It turns out that the singularities in 
the functions X^(fc) are intimately linked to the bound state spectrum and they 
ensure that the total number of fermion states remains the same. The issue is 
covered in detail in section 8, where a specific model is used to demonstrate these 
points explicitly. 

Consider first a Type 1 singularity. Denoting the the last solution of the 
discretisation condition before the singularity by and the first solution above 
the singularity by in the case of a Type 1 singularity we see (fig.[|) that, 

cot-i(-X(A;,^J) |p~7r, (157) 

and, 

cot-i(-X(fci^J) |p~ 0. (158) 
Using is+ = is- + 1 and the discretisation condition (|154[), we see that. 



(fci,+ - V)^o = 0(-i-)+/3(z,^j7r-/3(v)7r. (159) 

As Ro becomes large we require that the momentum shifts are continuous, i.e. 
(fci^+ — ki^- ) ~ Thus we have a jump in /3; 

/3(z,+ )=/?(z,-) + l. (160) 

This is precisely the jump expected due to a Type 1 singularity skipping a branch 
of the cot function. The singularity increases the value of 5 by one, thus we expect 
(3 to increase by one. This jump ensures continuity of the discretisation condition 
about a Type 1 singularity. If we start at large momentum and work towards 
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small momentum, on passing the Type 1 singularity the inverse cotangent jumps 
upwards by vr, but at exactly the same point the function (3{i) jumps down by tt: 
hence there is no overall jump in the function and continuity across the singularity 
is ensured. 

A similar analysis for a Type 2 singularity gives, 

(3{ts+)=f3{ts-)-l. (161) 

Again this jump is exactly that expected from a Type 2 singularity giving two 
solutions on one branch of the cot function, it ensures continuity of the discreti- 
sation condition around a Type 2 singularity. 



It can be seen from ( 160 ) and ( 161 ) that the physical requirement of continuity 
in the spectrum arises naturally in this formulation. 

Having fixed the discretisation condition uniquely for all values of the mo- 
mentum the next step is to use the discretisation condition to formulate the reg- 
ularised trace over the continuum states defined by (|50D. The manner in which 
the discretisation condition has been set up allows an integral over the continuum 
states to be defined in terms of the scattering data, found by solving the Dirac 
equation. 

The discretisation condition can be written uniquely as, 

kiRo = k^Ro + A^{ki), (162) 
where the free momentum is defined by, 

k^Ro = ivr, 

and, 

A^(fc) =p^ + e, + cot-\-X^{k)) I™"*, (163) 

where i G {1, 2, } and cot~^(— X^(A;)) I™"* is continuous and takes the princi- 
pal value for k larger than the last singularity in X^{k). It should be noted that 
for the pure Z-string there are two discretisation conditions of this type for each 
angular momentum mode labelled by ± in ( p.63| ). These arise from the two roots 



of the quadratic equation (|147|) which had to be solved in the process of setting 



up the discretisation condition for the Z-string background. 



Appendix B: Computation of the Z-string Heat Kernel Coefficients 

The purpose of this appendix is to provide details of the calculation of the 
coincidence limit heat kernel coefficients required for the computation of the 
counter term in the pure Z-string background. 

Tr[ai]. Using the standard recurrence relations we have, 

Tr[ai] = -a(x) = -4g^f{r)-2iq{aL-aR)da,+iq{aL-aR)diZ, TraV\ (164) 
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Using the standard anti-commutation relation satisfied by the Pauh matrices, 
{a^,a^} = 26^\ gives, Tra^a^ = 25'-', hence the second and third terms in ( |164| ) 
cancel, leaving, 

Tr[ai] = -^'f{r). (165) 

At this juncture it is however useful to note that individually the second and 
third terms in ( p.64| ) vanish because the divergence of the gauge field in the string 
background is zero, that is diZi = 0, which is easily shown by direct calculation 
using (0). 

Tr[a2]. Again using the standard recurrence relations we find, 

Tr[a2] = ^rrr,,r,, + iW(x). (166) 

The first term in ( |166| ) is straightforward to compute: Tij is the standard field 
strength tensor; 

Tij = diTj-djT, + [Ti,T,]. (167) 



In the case of the pure Z-string the third term in ( [L67| ) is zero because of the 
diagonal nature of Fj, hence, 

r,, = a.r, - a,r, = _° ^ j , (les) 

where, 

% = diZj - djZ,. (169) 

This is expected because the Z-string is really a Nielsen-Olesen U{1) string em- 
bedded in the U{l)z subgroup generated by the Z gauge bosons. Clearly Zij is 
anti-symmetric and so, 

r,T., = -2g^Z,^, ( f ^O^)- (170) 

Using the specific forms of the gauge fields in the string background (|^) implies 
that, 

= (171) 
r 

and therefore. 



z/'(r)' 



2 



Trr,,r,, = -4(a^ + ai)M:^ . (172) 



r 



As we only need its trace, only the diagonal elements of a(x)^ need to be 
considered. Using diZi = 0, it can be shown that, 

TraixY = 4g^f{rY-q\al + al)TT{a^a%Z,y 
+2g%{f{r)e^^')d,{f{r)e^^') 
+2q^g^f{rfZ,Z, 

±2zqg'f{r)Z, (e±'^9.(/(r)e^'^) - e^^'d,{f{r)e^'')) . (173) 
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Now the explicit forms of the gauge fields given in (0) can be substituted into 
(|173|) . The computation is straightforward but it is useful to give the explicit 
forms of some of the individual terms which appear in (IT] 



' v{r) 



qr 



.lAr) 
qr 



e-*' 
e'^ 



e"^' 
-e'^ 



'fir) 



(174) 

(175) 
, (176) 
(177) 



Using (|174|) to ( |177| ) in (|173|) gives the required expression for Tra(x)^ expressed 
purely in terms of the Nielsen- Olesen profile functions and the parameters in the 
theory. After a little algebra it can be shown that, 



+2/ (fir))' + 



'fir) 



+2g 



2fir)Mr)' , , 2fir)Mr) 



(178) 



It is then easy to show by direct substitution of ( |172| ) and ( p,78|) into (|166| ) that 
the expression for Tr[ai] is indeed that given by (|6^). 
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